Magneto-optical properties of multilayer graphenes 
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The magneto-optical absorption properties of graphene multilayers are theoretically studied. It is 
shown that the spectrum can be decomposed into sub-components effectively identical to the mono- 
layer or bilayer graphene, allowing us to understand the spectrum systematically as a function of the 
layer number. Odd-layered graphenes always exhibit absorption peaks which shifts in proportion to 
VB, with B being the magnetic field, due to the existence of an effective monolayer-like subband. 
We propose a possibility of observing the monolayer-like spectrum even in a mixture of multilayer 
graphene films with various layers numbers. 



I. INTRODUCTION 

The unusual electronic property of the atomically thin 
graphene films has been of great interest. Recently the 
optical absorption spectra were measured in graphene- 
related systems under magnetic fields f^i^i^i^i^i^ In this pa- 
per we theoretically study magneto-optical spectra of the 
graphene multilayer. 

The monolayer graphene is a zero-gap semiconduc- 
tor with the linear dispersion analogous to the zero- 
mass relativistic particle. In presence of a magnetic 
field it gives an unusual sequence of the Lan- 
dau levels with spacing proportional to ^fB in both 
of the electron and hole sides. The transport prop- 
erties in such a unique band structure were studied 
and found to be significantly different from the conven- 
tional systemj^i^iiSiiiiii^ii^ Recent experimental realiza- 
tions of monocrystalline graphene opened the way for di- 
rect probing of those unusual properties . ^"^'^^i^^ The Lan- 
dau level structure of the single-atomic sheet of graphene 
was investigated through the quantum Hall effec t^^i^^ 
and the cyclotron resonancei^ The optical response in 
graphene monolayer was theoretically studiediiiiiiii^ 

The multilayer systems containing few layers of 
graphene are also fabricated,— have attracted broad 
interest as wellf^S, There the interlayer coupling drasti- 
cally changes the structure around the band touching 
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On the other hand, recent 



observations of the magneto-absorption spectra of thin 
epitaxial graphite^ show i/B-dependent transition peaks 
just as in monolayer graphene. ^'^ Similar evidences for 
the linear dispersion were also found in thicker graphite 
systems Quite recently the cyclotron resonance 
was measured in graphene bilayer.^ In theories, the elec- 
tronic structure in magnetic fields has been extensively 
studied for three-dimensional (3D) graphite^'^'^i^'^'^ 
and for few-layers graphenes .^^i^"^ The optical ab- 
sorption was theoretically investigated for the bilayer 
graphene i^i^i 

Here we study the optical absorption properties of the 
AB-stacked multilayer graphenes in magnetic fields sys- 
tematically as a function of the layer number. We decom- 
pose the Hamiltonian into subsystems effectively iden- 



tical to monolayer or bilayer graphene )^ and express 
the spectrum as a summation over each of them. We 
present in Sec.|II]the Hamiltonian decomposition and the 
Landau-level structure of the multilayer graphene as well 
as the formulation of the optical absorption. We show 
the numerical results in Sec. IIIII and discussion in Sec. 
llYl 



II. FORMULATION 

We consider a multilayer graphene composed of N lay- 
ers of a carbon hexagonal network, which are arranged 
in the AB (Bernal) stacking. The system can be de- 
scribed by a k-p Hamiltonian based on 3D graphite 
model.— i^i^ The effective models were derived for the 
monolayer grapheneji'^'^^i^ the bilayerfSi and the tri- 
layer and morcj ^^'^^ For the simplicity we include the 
nearest- neighbor intra-layer coupling parameter 70, and 
the inter-layer coupling 71 between A and B atoms lo- 
cated vertically with respect to the layer plane. The 
band parameters were experimentally estimated in bulk 
graphite as 70 ~ 3.16 eV^S and 71 « 0.39 eV.^ The 
effects of other band parameters neglected here will be 
discussed in Sec. IIVI 

The low energy spectrum is given by the states in the 
vicinity of K and K' points in the Brillouin zone. Let 
\Aj) and \Bj) be the Bloch functions at the K point, 
corresponding to the A and B sublattices, respectively, 
of layer j. For convenience we divide carbon atoms into 
two groups as 



Group I 
Group II 



Bi, A2, B3, 
Ai, B2, A^, 



(1) 
(2) 



The atoms of group I are arranged along vertical columns 
normal to the layer plane, while those in group II are 
above or below the center of hexagons in the neighboring 
layers. The lattice constant within a layer is given by 
a = 0.246 nm and the distance between adjacent layers 
Co/2 = 0.334 nm. 

If the basis is taken as \Ai),\Bi); \A2),\B2); 
l^jv), \Bn), the Hamiltonian for the multilayer graphene 
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around the K point becomes 



H = 



/Ho V 

yt Ho Ft 

V Ho V 

\ 



(3) 



with 



Ho 





VTTa 



V7T- 





V 



■J 




71 



(4) 



where 7r± = tt^; ± itTj, with tt 



-i/iV + eA, the vec- 



tor potential A, and v is the band velocity of monolayer 
graphene, which is related to the band parameter via 
V — \/3a'-fo/2h. The effective Hamiltonian for K' is ob- 
tained by exchanging 7r_|_ and 7r_ . 

The Hamiltonian ([3]) can be decomposed into smaller 
subsystems for the basis appropriately chosen. First, 
we define the orthonormal sets 



|0™) = V/(i)|i?i> 
|0|"))=V^,(i)|Ai) 

where 



with 



iV+ 1 



sinjK;, Ki 



>i(3)|^3> 



Itt 



2(iV+ 1)' 



1 = ~{N~1), ~{N~3),...,N~l. 



(5) 



(6) 



(7) 



Here, I is an odd integer when the layer number N is 
even, while I is even when N is odd, and therefore I — 
is allowed only for odd N. 

Next, for to > 0, we take the basis 



{(I0^"^> + l0^"l»/^/2, (I^W) + \<f>^ll))/V2, 



(8) 



For TO = 0, we take the basis {|0o"^), |0o^)}. Then, 
the Hamiltonian has no off-diagonal elements between 
different to's. For m > 0, the sub-Hamiltonian within 
the basis Eq. ^ becomes 



with 



vn^ 

A„.7i 

A„i7l VTT- 

VTr+ 



2 cos Kr, 



(9) 



(10) 



which is equivalent to the Hamiltonian of bilayer 
graphene, while the inter-layer coupling 71 is multiplied 
by Am- For to = 0, we have 



U7r_ 
m+ 



(11) 



which is identical to the Hamiltonian of the monolayer 
graphene. These subsystems are labeled as 



TO 
TO 



:0,2,4, 

1,3,5, 



■ ,N- 



(odd N), 
(even N). 



(12) 



The eigenstate of a finite-layered graphene can be re- 
garded as a part of a standing wave in 3D limit, which 
is a superposition of opposite traveling waves with ifc^. 
The quantity k {— Km) in our representation corresponds 
to the 3D wave number via k — \kz\co/2. Thus the 
monolayer-type subband k = 7r/2 is related to a H point 
in the 3D Brillouin zone, while no states exactly corre- 
spond to ^2 = since k never becomes zero. 

The Landau levels of the monolayer-type states are 
given by 



(13) 



with n = 0, 1, . . . and s = ±, where s = + and — repre- 
sent the electron and hole bands, respectively, and only 
s = -f is allowed for n = OJ Here As is the magnetic 
energy, defined by 



Ab ^ V2hv^eB. 



(14) 



The Landau-level structure of the bilayer-type Hamil- 
tonian ([9]) was obtained previously'*® and can be analyt- 
ically derived by noting that 7r± are associated with the 
ascending / descending operators of the Landau levels^! 
in a similar way to that for 3D graphitei^!^^ The eigen- 
function can be written as 



(15) 



with n > ~1 and amplitudes Ci. Here ipn,kix,y) is the 
wavefunction of the nth Landau level in conventional 
two-dimensional system, given in the Landau gauge A = 
(0,^2:) by ifn^k = i"(2"n!v^/)-i/2e*'=2'e-^'/2iJ„(z) with 
z ^ {x + kP)/l and being the Hcrmite polynomial. 
We define ipn,k = for n < 0. 

For n > 1, the Hamiltonian matrix for the vector 
(ci, C2, C3, C4) then becomes 



/ AbV^ 

AsVn A71 

A71 
\ OA 



AB\/n+ 1 





(16) 



where the index of A„i is dropped. This immediately 
gives four eigen values 



V2 



(A7i)2 + (2n+l)A 



-f/iV(A7i)4 + 2(2n + l)(A7i)2A| + A* 



1/2 



, (17) 



where /i = ± correspond to the higher and lower sub- 
bands in the limit of zero magnetic field, respectivelyi^ 
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In the following we use the notation jj. = H,L instead of 
+, — to avoid the confusion with s — ±. The eigen states 
can be labeled by n, /z, s, and k. 

For n = 0, the first component of the wave function 
T5l) disappears and we have only three levels, 



0, 



7? 



(18) 
(19) 



At n = — 1 only the last component survives in (|15p so 
that we have only a single level in the lower subband, 
E-i^L = (the level e_i,_f/ does not exist). 

In small magnetic fields, the Landau levels for the lower 
subband in the region e <C A71 ar e approximately given 
by Sn,L,s ~ s{heB/m*)y/n(n + 1) with the effective mass 
m* = A7i/(2v^).^i Thus the level spacing shrinks much 
faster in i? ^ than that in the monolayer oc \/B. 
The ratio of the first gap of the bilayer-type subband, 
heB/m*, to that of the monolayer, A^, is given by 
Ab/(A7i). 

Figure [T] shows the Landau levels of the bilayer-type 
Hamiltonian as a function of k, in the magnetic field given 
by Ab/ji — 0.5. The bilayer levels become those of two 
independent monolayers at k — tt/2, where the effective 
inter-layer coupling A71 vanishes. The levels become flat 
around k — 0, where dX/dn vanishes. In the bottom 
panel we show the list of k for every layer number iV. 
The top and bottom panels share the horizontal axis; the 
Landau levels in the specific point in the bottom panel 
are shown directly above. 

The velocity operator for the sub-Hamiltonian Tim is 
given by — ~{i/h)[x,Hm] = dHm/dnx- There are no 
matrix elements connecting different m's. For bilayer- 
type subband, Vx has a non-zero matrix element only 
between the Landau levels with n and n±l for arbitrary 
combinations of fi ~ L and s = ±. This is explicitly 
written as 

{n',fi',s';k'\vx\n,^i,s;k) ^ vSkM' [ic'^C2 + c'^C4)S„^n'-i 

+ {c'2*Ci+c'^*C3)Sr,,n' + l] , (20) 

where Ci and are the eigenvectors of the matrix 
(fTB)) . corresponding to the Landau levels {n,^,s) and 
{n' , fj,\ s'), respectively. For the monolayer-type band, 
we have 

{n',s';k'\vx\n,s;k) = vSk^k' [ci*C2^„,„'-i + C2Ci6n,n'+i] , 

(21) 

where (ci, C2) is (0, 1) for n = and (s, l)/\/2 for n > 1 
for if -point 4 

To estimate the optical absorption intensity, we calcu- 
late the real part of the dynamical conductivity a^i-xiuj). 
The relative transmission of the sheet to the vacuum, for 
the linearly polarized light incident perpendicular to the 
plane, is related to this quantity viai2^ 
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FIG. 1; (top) Landau levels of the bilayer-type subband as 
a function of k with A = 2 cos k. Magnetic field strength 
is taken as As/71 = 0.5. (bottom) Lists of k in Af- layered 
graphene. Empty and filled circles represent even and odd 
TV's, respectively. 



As will be shown below, the expansion is valid except 
in thick multilayer graphenes for which the absorption is 
significant. The dynamical conductivity can be written 
in usual manner as 



a, (3 



6/3 + hiu + iS 



,(23) 



Rea^xito). (22) 



where S is the area of the system, is the velocity op- 
erator, 6 is the positive infinitesimal, /(e) is the Fermi 
distribution function, and \a) and describe the eigen- 
state and the eigen energy of the system. 

In the simplest approximation, we include the disorder 
effect by replacing d with the phenomenological constant 
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h/r = 2T and taking the ideal eigenstates as a, /3. The 
conductivity can then be written as a summation over 
all the contributions of the subsystems, which are inde- 
pendently calculated. Correspondingly, we compute the 
density of states per unit area as 



D{e) = -^Im 



1 



(24) 



with the ideal eigenstates a. 

The dynamical conductivity at zero magnetic field 
was calculated for the monolaye r^^'^^ and the bilayer 
graphene.— For the ideal monolayer at ep — 0, the 
expression apart from uj — Q becomes a frequency- 
independent valueiiiii 



9v9s e' 
16 h' 



(25) 



with g„ ~ 2 being the valley (if, K') degeneracy and 
cjs — 2 the spin degeneracy. Note that the dynamical 
conductivity has a singularity at {sp, uj) = (0, 0), which is 
removed if level-broadening effect is included properly. 
The expression for the effective bilayer Hamiltonian ([9]) 
with Ei? = is given bySS 



ReCTj,j;(cj) = 

huj — 2A71 
hw — A71 



9v9s e 
16 h 



e{hLU - A71) 



e{huj - 2A71) 



huj + 2A71 



huj + A71 



(26) 



where d{x) = 1 for x > and d{x) = for a; < 0. The 
above shows that the typical value of the real part of the 
conductivity for fiio/^i < 1 is {gvg s/ 16) {e^ /h) per layer. 
By noting that /he k, 1/137, we see that the expansion 
in Eq. ^ is valid roughly for < 20. 



III. NUMERICAL RESULTS 

Figure[2]shows the plots of Re a^x {'^) for the monolayer 
and bilayer graphenes in several magnetic fields. Here we 
take r/71 = 0.01, ep = 0, and zero temperature. Dotted 
lines penetrating panels represent the transition energies 
between several specific Landau levels as a continuous 
function of . The peak positions of each panel corre- 
spond to the intersections of those and the bottom line 
of the panel. 

In the monolayer the peak position obviously shifts in 
proportion to ^/B (i.e., cx As). In the limit of vanishing 
magnetic field, the conductivity eventually becomes the 
value given by Eq. (^5]) . The spectrum in the bilayer is 
rather complicated; starting from w = 0, we first see the 
series of the transition peaks within L bands from lj — Q, 
and then those between L and H enter for hu) > 71 and 
lastly those within H bands for hw > 2ji . The every 
peak position behaves as a linear function of B (cx A^) 
in weak fields but it switches over to -s/S-dependence as 
the corresponding energy is going out of the parabolic 



band region. In small fields the peaks are smeared out 
more easily in the bilayer than in the monolayer. The 
conductivity converges to the zero-field curve with a step- 
like structure at e = 71, which is expressed as Eq. ([^^ 
in the clean limit. 

Figure [3] shows the plots of Kea-xxi^^) from iV = 1 to 5 
with two different magnetic fields with Ab/7i =0.1 and 
0.2. We again take r/71 — 0.01, ep = 0, and zero tem- 
perature. The results are shown separately for each sub- 
band. In every odd layers the monolayer-type subband 
gives the identical spectrum. All other bilayer-types give 
different spectra depending on A. The quantized feature 
is more easily resolved in a subband with a smaller A, 
because of its narrower level spacings. In zero field limit, 
every bilayer- type spectrum has a step at £ ~ A71 , where 
the excitation between L and H bands starts. 

It is intriguing to consider how the absorption spec- 
trum looks like when the sample is a mixture of thin 
graphene films with various layer numbers. One might 
think the discreteness of k is easily smeared out and 
we just get the 3D limit, but it is not always the case 
as we will show in the following. We here calculate 
the dynamical conductivity averaged over the samples 
TV = 1,2,3, • • • ,20. We show in Fig.Hplots oiReaxx{i^) 
for different magnetic fields with Ab/7i = 0.1, 0.2, and 
0.3 and in Fig. [Sla gray-scale plot of Heaxxi^, ^b)- 

Surprisingly we still see the series of peaks in the 
monolayer graphene cx ^/B. This comes from the 
monolayer-type subbands which appears in every odd 
layered graphene. The visibility of the monolayer-type 
signal depends on the ratio of the number of monolayer- 
type subbands to the total; in the present case, this is 
10 to 110. We expect that the signal of monolayer grad- 
ually becomes invisible as the maximum layer number 
-^max becomes larger, because the total subband number 
increases as cx N^^^ while the number of monolayer-type 

as cx A^'max- 

We have another set of dominant peaks, which can 
be identified as the bilayer- type with k = (A = 2). 
Although there is no subband which exactly takes this 
value, many subbands around k ~ have almost the 
same peak positions as the Landau level is flat against k 
there and gives similar spectra. Every peak shifts upward 
with respect to the original position of k = 0, since the 
Landau level spacing is generally wider for larger k. Un- 
like the monolayer-type signal, this would survive even in 
the 3D limit, since the finite region in k (not a point) can 
contribute to this spectrum. When decreasing the mag- 
netic field, however, the bilayer-type peaks are immedi- 
ately blurred due to rapid _B-linear dependence, while the 
monolayer peaks survive even in relatively smaller mag- 
netic field. In zero-field limit, we are left with a bump at 
huj — 271, which comes from the H-L transition step of 
the bilayer-type subbands with k ~ 0. 

Just in the same way as the monolayer-type subband 
(k — n/2) appears in every two layers, the bilayer-type 
subband with k = 7r/3 enters in every three layers (iV = 
2,5,8, •• •) and that with k = 7r/4 in every four layers 
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FIG. 2: Real part of the dynamical conductivity of the mono- 
layer (top) and bilayer (bottom) graphenes plotted against the 
frequency ui, calculated for different magnetic fields (specified 
by A_g) and r/71 = 0.01. Dashed curves indicate the transi- 
tion energies between several Landau levels in the ideal limit. 
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(iV = 3, 7, 11, • • • ). We can see the H-L transition peaks 
of those k's in Fig. [5l while L-L peaks are hidden by 
other dominant contributions. 

The similar analysis is available for the density of states 
(DOS). In Fig. [HI the top panel shows DOS averaged over 
iV = 1, 2, 3, • • • , 20 as a function of the Fermi energy. The 
bottom panel shows the corresponding plot for the local 
density of states (LDOS) on the top layer, defined by 
the number of states per unit energy width and per unit 
area on the layer. We also present in Fig. [7] the two- 
dimensional plots of DOS and LDOS on (e, As)-plane, 
where the gray-scale shows the relative value from the 
zero magnetic field. 



FIG. 3: Real part of the dynamical conductivity of A'' = 
1,2, ■■■ , 5-layered graphenes as functions of the frequency, 
calculated for two different magnetic fields A_b/7i =0.1 and 
0.2. We take r/71 = 0.01. 



In DOS, we observe the several peaks coming from 
the monolayer-type subband similarly to the optical ab- 
sorption spectra. The peaks from the bilayer-type sub- 
band with K = become prominent in the high-field re- 
gion Ab/7i>0.2. In LDOS, interestingly, the peaks of 
the monolayer-type subband are much more pronounced. 
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FIG. 4: Plots of Re <Jxx{io) averaged over the layer numbers 
from A'' = 1 to 20, for magnetic fields with Ab/^i = 0.1, 0.2 
and 0.3. Vertical solid and dashed lines represent the ideal 
transition energies for the monolayer-type (k — it/2) and the 
bilayer-type {k = 0) subbands, respectively. 



while those of k = are strongly suppressed. This can be 
understood by the wave function defined by Eq. ([H]). If we 
look at a state in the subband m in A^-layered graphene, 
the wave amplitude on the top layer {j = 1) always ac- 
quires the factor 7/'m(l) = \/2/(iV +1) sin Km- Obviously 
this takes maximum in the monolayer-type (k — 7r/2) 
and zero at k = 0, and thus the monolayer-type state 
contributes the most to the surface LDOS. The bilayer- 
type signals of the subbands k = tt/3 and 7r/4 are also 
visible in LDOS while they are hidden by k = in DOS. 



IV. DISCUSSION 

Recently the optical absorption spectrum was mea- 
sured in the epitaxial thin graphite films and the 
monolayer-like signal was observed, while the detail 
profile of the system remains unclear. Similar ^/B- 
dependent features were also observed in the samples con- 
taining high-number of graphene layers 100) grown on 
SiC substrate,^ and in a thin graphite sample of thick- 
ness ~ 100 nm exfoliated from highly-oriented pyrolytic 
graphite.— Those results are nontrivial because, if the sys- 
tem is a real three-dimensional bulk graphite, the spec- 
trum would be contributed mainly from the states around 
fc^ = (k = in our discussion) where the Landau levels 




2 rico / Y, 



FIG. 5: (Top) Density plot of Re(Ti:i(a;, A_b) averaged over 
the layer numbers from A'^ = 1 to 20. F is set to O.OI71. (Bot- 
tom) Corresponding plot for the transition energies between 
the several Landau levels. 
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Q 3 0.5 



0.8 E / y, 



FIG. 6: (Top) Density of states averaged over the layer num- 
bers N — 1 to 20, at As/71 = 0.3. (Bottom) Corresponding 
plot for the local density of states on the top layer. Vertical 
solid and dashed lines indicate the ideal Landau level ener- 
gies for the monolayer-type {k — n/2) and the bilayer-type 
(k — 0) subbands, respectively. 



7 



DOS difference 
(units of g^gj^fijtfi'^v^) 




0.2 0.4 0.6 0.8 e/y, 



LDOS difference 
(units of g„gj,/27iti=^v=} 




0.2 0.4 0. 




FIG. 7: (Top) Density of states (measured from the zero- 
field value) averaged over the layer numbers from = I to 
20, plotted in (e, AB)-plane. T is set to O.OI71. (Middle) 
Similar plot for the local density of states on the top layer. 
(Bottom) Ideal Landau level energies corresponding to several 
dominant peaks in above two panels. 



are flat with respect to kz . One possible scenario for this 
is that the system can be regarded as a compound of mul- 
tilayer fragments with various small layer numbers, and 
the monolayer-like spectra of all the odd layers are ob- 
served. It should also be mentioned that the local density 
of states on the surface of graphite was observed in the 
experiment j ^^i^^ Our calculation predicts that the pro- 
nounced monolayer-type spectrum would be observable 



in a multilayer graphene. 

While we adopted a simplified effective-mass model 
in which only 70 and 71 are included, here we briefly 
mention the effects of other hopping parameters. The 
parameter 73 neglected here couples group II atoms on 
neighboring layers. This is responsible for the trigonal 
warping of the band dispersion, but gives only a slight 
shift in the Landau level energies except for the low 
energy region (< 10 meV)M^^ Therefore, it would 
hardly affect the peak positions in the absorption spec- 
tra while may modify the amplitudes through the ma- 
trix element changes. The parameter 74 couples group I 
and II atoms sitting on the neighboring layers, such as 
Aj ^j+i or Bj Bj^i. This parameter introduces a 
small electron-hole asymmetry in the band structure, but 
does not change the qualitative feature of the low-energy 
spectrumi^^ 

We also neglected the vertical hopping between the 
second-nearest neighboring layers for group II and I 
atoms, which are parameterized by 72 and 75, respec- 
tively. Including those parameters mainly shifts the 
zero energy (the band touching point) upward or down- 
ward, depending on each subsystemi^i2£i22 in the tight- 
binding model^i21 and the density functional theory^- 
estimate the shift SE at the order of 10 meV. In 3D 
limit, this corresponds to the band dispersion along fc^- 
directionii2ii2^i^i2ii2^ The zero-energy shift leads to the 
electron or hole doping, and gives a change of the absorp- 
tion spectrum in the region Huj < 2SE. 

The effective mass model is no longer valid when the 
energy is as high as the intra-layer coupling 70 '--^ 3 eV. 
The lattice effect appears as trigonal warping in the band 
dispersion in higher energies, ^^'^"''^^ while this should be 
distinguished from the trigonal warping discussed above, 
which is due to the extra band parameter within the ef- 
fective mass model. The frequency region covered in our 
calculation, h(jj<2.5ji ~ 1 eV, roughly corresponds to 
the energy region |e|<0.5 eV. The deviation in eigen en- 
ergy is estimated at 5% at e = 0.5 eV and can be treated 
perturbationally,^^ although it grows as the energy in- 
creases out of this region. This anisotropy constitutes a 
major part of the chirality dependence of optical spec- 
tra in carbon nanotubes, enabling the assignment of the 
structure of individual nanotubes. 

Lastly, while our model is based on the bulk 3D 
graphite, it should be noted that the band parameters 
in few-layered graphenes are not exactly the same as 
those for the bulk graphite, but generally vary depend- 
ing on the layer number. There is a theoretical at- 
tempt to obtain accurate electronic structures for few- 
layered graphenes, using the density functional theory 
with the local density approximatioui^^ The calculation 
beyond the local density approximation was also pro- 
posed, which properly treats nonlocal van der Waals in- 
teraction coupling graphene layers in the density func- 
tional framework)^ The study of the optical absorption 
in a refined band model is left for a future work. 

In conclusion, we have presented a systematic study 
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of the optical absorption properties and the density of 
states in the multilayer graphenes as a function of layer 
numbers. The spectrum can be understood through 
the decomposition into sub-components, each of which 
is equivalent to the monolayer graphene or the bilayer 
graphene with single parameter k. We proposed that the 
monolayer-like spectra is possibly observed in the mix- 
ture of the multi-layered graphene, contributed by the ef- 
fective monolayer subbands existing in every odd-layered 
graphene. 
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